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1.  Introduction 

Our  understanding  of  nonlinear  dynamics  of  deep  water  gravity  waves 
has  grown  substantially  in  recent  years.  We  feel  that  the  lion's  share  of 
this  progress  should  be  attributed  to  the  staff  of  the  TRW  Fluid  Mechanics 
Department.  Most  of  their  findings  are  summarized  in  an  extensive  review 
article  by  Yuen  &  Lake  (1982),  which  served  as  our  main  reference.  Much  of 
this  progress  is  based  on  applications  of  the  so-called  Zakharov  equation 
which  was  originally  derived  by  Zakharov  (1968)  for  infinitely  deep  water. 

In  order  to  extend  the  range  of  application  we  re-derive  Zakharov's  equation 
for  finite  water  depth  (in  section  2)  and  show  its  relations  to  the  cubic 
Schrodinger  equation  and  to  Hasselmann's  nonlinear  interaction  model  (in 
section  3).  It  is  generally  accepted  now  that  the  Zakharov  equation  is 
superior  to  all  other  existing  approximate  models  as  far  as  Class  I  inter¬ 
actions  are  concerned. 

The  term  'Class  I  interactions'  refers  to  nonlinear  interaction 
processes  at  the  lowest  possible  order;  for  surface  gravity  waves  this 
occurs  at  third  order  in  the  nonlinearity  parameter  e.  Generally  speaking. 
Class  I  Interactions  require  the  coexistence  of  resonating,  or  nearly 
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resonating,  wave  quartets.  The  time  scale  of  Class  I  interactions  is  e  P 
where  P  is  a  typical  wave  period. 

The  structure  of  the  surface  gravity  wave  dispersion  relation  does  not 
enable  nonlinear  interaction  at  shorter  time  scales  (e  P)  which  occur  in 
many  other  physical  systems, (e. g. ,  capillary  waves). 

While  Class  I  interactions  are  basically  four  wave  interactions,  the 
special  case  where  one  of  the  waves  is  taken  into  account  twice  so  that 
only  three  waves  are  considered,  has  attracted  much  attention.  These  cases 
which  lead  to  what  sometimes  is  called  Benjamin-Feir  instabilities,  display 
many  of  the  features  of  the  more  general  quartet  interaction.  Interactions 
including  a  smaller  number  of  waves  —  as  two  waves  each  taken  into  account 
twice,  or  one  wave  taken  into  account  four  times  —  are  also  possible,  but 
display  a  degenerated  type  of  interaction  which  manifests  itself  in  Stokes- 
type  second  order  corrections  of  the  frequency  (see  Longuet-Higgins  and 


Phillips  1962).  Numerical  linear  stability  analysis  of  the  exact  finite 
amplitude  Stokes  wave,  by  McLean  (1982a,b),  as  well  as  experimental  evidence 
by  Su  et.  al .  (1982)  and  Su  (1982),  reveal  the  importance  of  Class  II  inter¬ 
actions,  which  are  basically  quintet  interactions.  These,  much  less  studied 
interactions,  occur  at  fourth  order  in  e  and  have  a  typical  time  scale  of 
e  3P.  Nevertheless,  for  high  enough  steepnesses  McLean's  study  shows  that 
Class  II  instabilities  become  dominant.  Here  again,  three  waves  —  one  of 
them  is  taken  into  account  three  times  —  form  a  nearly  resonating  quintet 
and  display  many  interesting  features.  In  the  second  half  of  section  2  we 
extend  the  derivation  to  fourth  order  and  derive  a  modified  form  of  the 
Zakharov  equation  which  accounts  for  both  Class  I,  and  the  higher  order. 
Class  II,  interactions. 

In  section  4,  we  use  this  equation  to  study  the  linear  stability  of  a 
uniform  w*ve  train.  The  solution  of  certain  long  line  evolution  problems 
is  under  way  and  will  be  reported  at  a  later  stage. 


2.  The  governing  equations 


The  equations  governing  the  irrotational  flow  of  an  incompressible, 
inviscid  fluid  with  a  free  surface  are: 


*  0  ,  -h<Z  £n(x.,t) , 


nt+  (•?♦)  •  (VrO-^'O 

♦  t+  -j(?<}>)Z+  S*®0 
'  4>z  «  r0  , 


2  *  n(x.»t) 


z  -  -h  . , 


(2.1) 


(2.2a,b) 


(2.3) 


where  0  is  the  velocity  potential,  n  Is  the  free  surface,  and  g  is  the 
gravitational  acceleration.  The  horizontal  coordinates  are  (xj,x2)  *  x_,  the 
vertical  coordinate  z  is  pointing  upwards,  h  is  the  mean  water  depth,  and 
t  Is  the  time. 

The- free  surface  boundary  conditions,  Eq.  (2.2), are  rewritten  in  terms 

of  $s  and  ws  ■  (•|^)|„»  the  velocity  potential  and  the  vertical  velocity 
pz  *  z*n 

component  at  the  free  surface,  respectively 


nt  +  (Vs)*(vxn)  -  ws[l  +  (vxn)23  »  o 


(2.4a) 


♦t  +  gn  +  j(vx$s)a"  ^  Cl+(vxn)l]«  o 


(2.4b) 


The  horizontal  Fourier  transform  of  these  equations  yields 
^(k  ,t)  -  ±  ||  (k.1-Jc2)$S(J^,t)  n^tjfi^-kj-k^d^d^-  ws  + 

CO 

*i],  ||j  (i2-k3)  ws(k.,t)nCk^t) n  (k.t)  e(k-krk2-k3)dk,dk2dk3.  0 


(2.5a) 


♦t  (jc,t}+gn(Jc ,t)-  —  jj  (k1-k2)«3(k1,t)<J>  (k2,t)6(k-kl-k2)dkidk2- 

— oo 
00 

“  -4—  J|  ws(k,  ,t)  ws(k2 ,t)  6  (k-k  1  -k2 )  dkj  d k 2  + 

-oa 

CO 

\\\\  (-iV"S{- I»t}  ws(k2,t)n(k3,t)n(kij,t)6(k-ki-k2-k3-j^) 


dk  dk  dk  dk,  r 
—l  -2  —3  -4 


(2.5b) 


where  the  two  dimensional  Fourier  transform  of  a  function  f(„x)  is  given  by 


?<J0  -  iff  f(x).-'S»d*. 


and  Dirac  6  function  is  defined  as 


*<*>  *  dr  l  *  -*2 


Taking  the  Fourier  transform  of  Laplace  Eq.  (2.1),  and  satisfying  the 
boundary  condition  at  the  bottom,  Eq.  (2.3),  gives 


$(k.»z,t)  «  $(k_,t)ch(|  k  ’  *+h)). 


(2.6) 


which  enables  one  to  write  $S  and  wS  in  terms  of  $U,t)  and  nU,t) 
as  follows 

00 

f  I  Ic  * 

4>S(x,t)»  ^  J  3(k.,t)[ch(  |  kjh)ch(  |J<|n(2L»t))+sh(  |kjh)sh(  |kjn(x,t))]  e  ~  *d£ 

~~  (2.7a) 

s  t“  *  ik*x 

w  (x>t)a  I  |k|4>(k_,t)[ch(|k|h)sh(|k_|n(x.,t))1+sh(|k^|h)ch(|Jt|n(x,t))]e  dk^ 

(2.7b) 
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The  next  step  in  the  derivation  is  to  express  wS  as  a  function  of  n  and  $s. 
This  is  the  first  step  which  requires  an  additional  physical  assumption. 
Assuming  that  jkjn  is  small,  we  pursue  the  following  procedure:'  (i)  replace 
sh(|kjn)  and  ch(jkjn),  in  Eqs  (2.7a,b),  by  their  Taylor  series  expansions 
up  to  order  (|k|n)3;  (ii)  express  n  by  means  of  its  Fourier  tranform  n  ; 
and  finally,  (iii)  take  the  Fourier  transform  of  Eqs.(2.7a,b) , 


♦  (k.t)«  $(Jc,t)ch(|kJh)+  JjlJ^lshClk^lhJo^.tJnCk^.tJfiCk,-^-^)  d^dk^  + 

+(2lip  Iff  V6"* li, lh>«<i, )dk|dk2d)c  + 


sh(  |k^h)»(kt  ,t)S(J^,t)n(l^.t)S(!cli.t)«(k.k|.Jtj^ -k^l 


•  dk  dk  dk  dk. 
“1  -2  “3  “** 


(2.8a) 


w 

S* (k.t)  -  |k|S(k,t)sh(|k|h)+  H-ljlkjl’chdkJhjI^.tJntJ^.tJStk-J^-J^JdJ^dk 


*sp  III  ¥ 


sh(|jL,jh)*(Jij  .tJnfJ^.tJn^.tJdCk-^-k^-k^JdJc^k^dl^  + 


(2^)3  ||f j  '"6^  ch(|k^Jh)$(k^  ,t)n(k^,t)n(k  ,t)Ti(k^,t)  (k-k^k^-k^-k^) 


dk  dk  dk  dk 
“1  ”2  “3 


(2.8b) 


Inverting  Eq.  (2.8a)  Iteratively,  in  order  to  obtain  5»  $  ($S),  *nd 
substituting  the  result  into  Eq.(2.8b),  yields 
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wS{k*,t)«  iJl) th(  |k.|h)$s(<k,t)  - 


'  27~  I!  l-l  I  th(  li!h)th(  li,  |h)-|jL,  |]$S  (k]  ,t)f(jc2,t)6(k-k.ik2)dk.1dk2- 


1~2  -1  -2 


-(SpJfJ  s(1^i*ii]*Js2.i3)$S(i1.t)n(k2,t)n{k3,t)6(k-k1-k2-k3)dkidk2dk  - 

-»«o  ^ 

CD 

'<*bIII s  (ii'ii,  ’iSj ^  *J^)  ^Jij  .  t)  q  ( ,  t)  n  ( ,  t)  n  ( ,  t )  6  U-kj  ) 


dk  dk  dk  dk,  (2  9) 

the  kernels,  S(1) (k.kj .k^) ,  S(2) (k.k,  ,k2,k  k^J ,  as  well  as  other 

kernels,  which  appear  throughout  the  derivation,  are  given  in  the  Appendix. 
Substituting  vr,  from  Eq.  (2.9)  into  Eqs.  (2.5a,b),  multiplying  Eq.  (2.5a) 

b»  r — a_ "I15  .  r  u(k)  I15 


|  2w(Jc )  |  ,  and  Eq.  (2.5b)  by  1  ["  “gg  -  | 

where  w(.k)  *  [g| Jc| th(  |k|h)],s  ,  (2.10) 

adding  these  equations  together,  and  defining  the  new  complex  variable 

b(k,t)  «  j-^T  ^  n  (k,t)  +  i  j~  to(Jl)  J*  $S'(£,t) 
yields  the  following  equation: 


bt  (Jc,*t)  +  iu(k)b(Jc,t)  + 


‘  lj[  V<1)(k,Jci,k2)b(kl,t)b(k2,t)6(k-ki^)dkjdJ^ 

-DO 

00 

iff  V(2)(k,k,.i2)‘'*(k|.‘»>(k2.t)6(!!+k)-k2)<ik|dk, 


1  Jj  V^,(k,kj,k2)b*(kj1t)b*(k2^)5(k+ki+k2)dkjdk2  + 

1  }j}  l)  dt,^  »t)b(k2,t)b(k^,t)6(jo-k>j-k^|^)dl^dki2dki^  h 

i  |j|  W(2> (k^k.,  .J^.jLjJb*^ *t)b(k2,t)b(k.3,t)6(^+k1-l^-^3)d^jdjc2djc3 
1  j|j  W(3>(Jk.lcJ ,k2,jc3)b*(k.l  .t:)b*(k2,t)b(k3,t)«(Jk4jLj4j^^)d^idJ^diL; 


1  jj]  W^4)(k,k1,k2,^)b*(kj,t)b*(j^,t)b*(l^,t)fi(k+ki+k2+J^)dk|dJ^d^  + 

1 1llK ])  (-*-i  »vvvb(-i  »t)b(^»tjb^*t)b^»tj6^,-vv^)d-idvvis4 + 

m 

1  III!  x(3)  *V!i3*Vbt^|  •t)b*(Vt)b(Vt)b(Vt)i(' VVV^t  Jdk^k^d^dk^ 

»  Ilf f  XW  (k.k,  .kj.kj.k^Jb^k,  ,t)b*( ,t)b*t ,t)b( <Bc , kB^dl^dk^ 

•O 

•  1  JJJf  X<^)  (k.»k.^  »1^|) b*(iL| »t)b*(k2»t)b*(j^,t)b*(l^,t)d(k>*Ht_^+k2+k^+kJf}dJkldk2dk^dk< 


Where  *  denotes  the  complex  conjugate 

The  relations  between  n»  4>S  and  the  complex  "amplitude  spectrum"  b  are 

i  o)(k),l/2 

n(k.t)  -  j^-|  [b(k,t)  +  b*(-k,t)] 

♦S(k,t)  -  “i  WiStItI  1/2  "  *>*C-k,t)3 


-  8  - 


We  assume  that  the  wave  field  can  be  divided  into  a  slowly  varying  (in  time) 
component  B  and  small  rapidly  varying  components  B',  B" ,  B"'  and  that  most 
of  the  energy  in  the  wave  field  is  contained  in  B.  These  assumptions  permit 
one  to  write 


b(k,t)  «  [e§U,t 2*3)  +  e2B‘  (k,t,t2,t3)  +  e3B"U,t,t2,t3)  + 


+  eV1  (k,t,t2,t3)]e‘i<1,(^)t  (2.13) 


where  c  is  a  small  parameter  representing  the  magnitude  of  nonlinearity,  and 

2  3 

the  slow  time  scales  are  defined  by  =  e  t,  t  =  e  t.  The  omission  of  the 
slow  time  t-j  ■  et  from  Eq.  (2.13)  results  from  the  fact  that  resonating  triads 
do  not  exist  for  surface  gravity  waves.  Substituting  b,  from  Eq.  (2.13),  into 
Eq.  (2.11)  and  arranging  the  terms  according  to  their  order  in  e  yields  the 
following  results: 


Order  c  -  is  satisfied  identically 


Order  e  -  gives  an  equation  for  B' 


181  .  f[  ,v(l)  s  s  .  e1(“-“1'“2.)t  +  v(2) 

9t  JJ  tV0,1,2BlB260-l-2  V0,1.2B1B2* 


i  (ur+a^-t^H 


0+1  -2 


+  V(3>  B*B*fi 
+  V0,l,2BlB260+l+2 


i  (w+u^+uOt 
e  L  dk^dkg  , 


(2.14a) 


where  we  have  introduced  the  compact  notation  in  which  the  arguments  k^.  in  V, 
B,  6,  a),  and  in  other  functions  in  the  sequel,  are  replaced  by  subscripts  i, 
with  the  subscript  zero  assigned  to  k^.  Integrating  Eq.  (2.14a)  with  respect 
to  t  and  keeping  *3  fixed,  gives 


B'  = 


•fKL 


|B, 


,2ur2u0-l-2  .  o>-u)1 

*  i  (u>+u>i-u>o)t 

(?)  -  c 

+  V0,1,2B1B260+1-; 


+  Vo!l,2BlB2S 


2  ur+u)-j  -(1)2 
i  (uj+u)-|  +^2^ 


0+1+2  UJ+Wn+l 


1+0)2 


}  dk^dkg  (2.14b] 


AA  ) 


Vr-.«V V 
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The  constant  of  integration,  which  corresponds  to  the  initial  phase,  has 
been  set  to  zero  without  loss  of  generality. 


Order  -  gives  the  following  equation 

1  1  IT*  Jj.  {^o]l,2,3  BlB2B360-l-2-3 


~(2)  ~ 

+  T0,l ,2,3  BlB2B36(Hl-2-3  e 


+  t<3)  i*ri  *  hkh^-ujU 

T0, 1.2,3  BlB2B3{0+1+2-3  e 


.  y(4) 

T0,l  ,2,3  BlB2B360+l+2+3  e 


i  ((j+^+u^+^H 


JdjOjdkg 

(2 


The  above  equation  consists  of  terms  of  two  types;  those  which  depend  on  the 
fast  time  t,  and  those  which  do  not.  This  enables  us  to  split  Eq.  (2.15) 
into  separate  equations: 


1  ^2*  JH  To!l»2.3  W3VI-2-: 


i  (urHj^-u^-WgJt 


djc-jd^dkg 


(2.16) 


3BM 

at 


,  ~n)  _  itw-UK-uu-toJt 

{T0,1,2,3  B1B2B360-.1-2-3  e  + 

"  ,;(2)  t(2)  ,  2*5  2  , 

'T0,l ,2,3  '  T0,l,2,3}  BlB2B360+l-2-3  c 

+  f(3)  6*5*5  ,  H— 

+  T0,l,2,3  BlB2B360+l+2-3  e 

*  ti.2.3  SX635W*2*3  e1(^^n,  ,2.17.) 


Here’  we  made  use  of  the  fact  that  the  only  exponent  of  e  which  may  become  zero, 
under  the  restriction  of  the  S  functions,  is  the  one  in  the  second  term  in 
the  r.h.s.  of  Eq.  (2.15).  This  fact  is  directly  related  to  the  definition 
of  a  nearly  resonating  quartet. 


6 0+1 -2-3  =  ^  *  | Wq+iu-j _w2_u3 1  1  0(e  ) 


(2.18a,b) 


Equation  (2.16)  is  the  so-called  Zakharov  equation,  with  the  kernel 


Tq  J  2  3’  ^or  near  resonance  quartets 


0,1 ,2,3  J  q  #  otherwise 


(2.19) 


used  as  a  mathematical  model  for  Class  1  nonlinear  interactions.  Integrating 
Eq.  (2.17a)  with  respect  to  t  gives  the  following  result  for  B", 

B”  =  'll)  {^ll,2,3  BlB2B360-l-2-3  6  + 

i  (0)+U).-UU-U)o)t 

+  (f(2)  _  t(2)  i  5*8  B  <5  e  1 

u0, 1,2,3  T0,1,2,3J  BlB2B360+l-2-3  w+u^-wg-wj 

i  (aj+u),+(jj0-u),)t 

-f3l  _  l  c  3 

+  o  „  ^ — — - + 

0,1, 2, 3  1  2  3  0+1 +2 -3  w+w-j+Wg-ajj 


~T4l 

+  T0,l  ,2,3  B1B2B36 


'  i  (uJ+W-j+OJg+UJ'jH 


0+1 +2+3  u+uu+i 


1+to2+u>3 


}  dtjdkgdkg 


(2.17b) 


Order  e 


1  1  TT-  -1  H^+  ||j|tiio!l),2.3,4'BrB2'Wo-l-2-3-4e 


i  (to— w-j  -u^-w^-w^)  t 


+  ^0  J,2,3,4BlB2B3B450+l-2-3-4e 


i  (ur+u-j-ujg-Wg-u^H 


+  Tj ( 3 )  ;*5*|  5  , 

U0,l,2,3,4BlB2B3B460+l+2-3-4e  + 

~t.)  ~*~*~*~  i(uH-w,+uu+cd,-a>.)t 

+  U0,l,2,3,4BlB2B3B46O+l+2+3-4e  + 


~f5l  '  -*~*~*~* 

+  U0,l,2,3,4BlB2B3B460+l+2+3+4e 


i  ( w+u^j  +052+u)2+w^  )  t 


}  dJli  djc^djc^djk^ 

(2.20) 
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In  order  to  split  Eq.  (2.20)  appropriately  into  two  separate  equations* 

dD  jvDll  1 

one  for  and  the  other  for  which  becomes  relevant  only  in  sextet 

interactions,  we  make  use  of  the  fact  that  only  the  second  and  third  integrands 
in  Eq.  (2.20)  enable  resonating  quintets.  Similarly  to  Eq.  (2.18),  the  nearly 
resonating  quintets  are  defined 


5 0+1 -2-3-4  *  0  *  InQ+ni"n2"n3“s4i  1  ) 


where  the  “Stokes  corrected"  frequencies,  are  given  by 


(2.21a,b) 


=  +  j  ejlTjljl  ^l'2  d-l  *  ejl  =  {1 


2  j  i  1 


j  - 1 


(2.21c) 


The  Stokes  corrected  frequencies  are  obtained  by  solving  Eq.  (2.16)  for 
degenerated  interactions,  namely:  "quartets"  formed  by  two  waves,  each  taken 
into  account  twice.  These  corrections  become  necessary  at  the  order  of  deriva¬ 
tion  considered  here.  Defining 

U0  1  2  3  4’  *or  near1y  resonating  quintets 


U, 


(2) 

0,1 ,2,3,4 


U 


(3) 

0,1 ,2,3,4 


0 

r.(  3) 


,  otherwise  (2.22a) 

[  °0°1  234’  for  near1y  resonating  quintets 

,  otherwise  (2.22b) 


we  obtain 


SB 

dU 


^  ^U0,l  ,2,3,4®!®2®3®450+l-2-3-4  e 


'  +  Uo!l,2,3,4glg2g36460+l+2-3-4  e 


i  (urh^-u^-u^-u^H 
i  (u+tu-i+a^-Wg-w^H 


)  dk-j  dkgdkjdk^ 

(2.23) 
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Finally,  the  two  orders,  Eqs.  (2.16)  and  (2.23),  are  combined' into  a  single 
equation  for  B  *  eB, 


if  -  III  T 


(2)  R*R  r  f  . 

0,1 ,2,3  BlB2B360+l-2-3  e 


i  (urh*>i  ui^Jt 


djc^d^gdk^  + 


HU  Uo!l  ,2,3,4  B1B2B3B4  6 

OD 

(3) 


i  (w+uV| 


0+1 -2-3-4 


dk^  dkgdkgdk^  + 


U 


Bt  30B0B  .6f 


i  (oj+uj^+oig-cj^-uj^Jt 


0,1 ,2 ,3 ,4  1  2340+1+2-3-4 


d^dkodMk, 


(2.24) 


Equation  (2.24)  is  a  modification  of  the  Zakharov  equation  that  accounts 
for  higher  order  interactions. 


3.  Comments  on  Zakharov  equation  for  finite  water  depth 
Denoting  e§  by  B  we  rewrite  Eq.  (2.16) 

3B(k,t) 


00 

fff  C>\  i  (w+u>,-u>.,-ui.3)t 

Jjj  ^i41»k2  43)B*(k1,t)B(J<2,t)B(k3,t)6(k+k.1-k2-k3)e  1  Z  3 

(3.1) 

The  first  order  free  surface  elevation  is  related  to  B  through  Eqs.  (2.12a), 
(2.13),  and  is  given  by 

i  7  fw(kh  1/2  ...  .. 

n(x»t)  *  2“  J  [-gg-j  {B(k.,t)el(- +  c.c}  dk  (3.2) 

•  00 

Equation  (3.1)  is  the  now  well-known  Zakharov  equation,  generalized  for 
water  of  any  constant  depth.  The  fact  that  Eq.  (3.1)  is  valid  for  finite  depth 
affects  only  the  expressions  for  u(k)  and  T^^k.k^  .kg,]^) »  which  become  depth 
dependent. 

The  purpose  of  this  section  is  to  show  the  connections  between  the  Zakharov 
equation  and  other  model  equations,  as  well  as  to  check  our  depth  dependent 
expression  for  T^.  Note  that  for  h  «•>  »  our  equation  for  T^,  Eq.  (A. 5c), 
gives  the  same  result  as  does  Appendix  A  of  Crawford,  et  al.  (1981).  This 
result  is  different  from  that  given  in  Yuen  &  Lake  (1982)  (even  after  correc¬ 
tions  of  minor  misprints).  This  apparent  discrepancy  is  related  to  the  special, 
almost-symmetric  (with  respect  to  kg  and  k^)  structure  of  the  Zakharov  equation. 
This  structure  allows  some  freedom  in  the  choice  of  T^,  i.e.,  T^ (k,^  ,ko,k,) 

- 1— t  — o 

can  be  replaced  by  aTv  '{k^k-j  ’kg’^)  +  ( 1  -a)  T'  jj^k^ »  with  arbitrary  a, 

without  altering  the  value  of  the  integral  on  the  r.h.s.  of  Eq.  (3.1).  Any 

T^,  obtained  in  some  legitimate  derivation  can  be  made  symmetric  in  k,,k-  by 

(21 

choosing  a  =  0.5.  This  symmetric  T'  ,  denoted  by  T,  is  a  uniquely  defined 
function  of  ]c,k^  ,kg ,k3  and  h  and  will  be  used  in  the  sequel. 

Relation  to  Hasselmann's  energy  transfer  model 


The  energy  transfer  equation  for  a  finite-depth  gravity-wave  spectrum, 
originally  obtained  in  Hasselmann  (1961),  was  rederived  by  Herterich  &  Hasselmann 
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(1980).  This  last  paper  served  as  a  reference  for  the  verification  of  the 
expression  for  T,  Eq.  (A. 5c).  By  reasonings  similar  to  those  in  Longuet- 
Higgins  (1976),  but  starting  from  the  Zakharov  equation  (instead  of  the  cubic 
Schrbdinger  equation,  used  by  Longuet-Higgins)  the  following  energy  transfer 
equation  is  obtained 

3C(k,t)  rrr 

- it~“  4ir  j  J  T(k,k1,k2,k3)[C(k2)C(k3)(C(k)  +  Cfe,))  -  C(k_)C(k1)(C(k2)+ 

-«o 

**■  C(Jc2))]6(Jc  -uij-tu^)  dk^dkgdk^  (3.3) 

O 

where  the  wave-action  spectrum  C  *  |B|  . 

For  strict  resonance  conditions,  which  are  implied  by  the  two  6  functions 
in  Eq.  (3.3),  T  is  also  symmetric  in  its  two  first  arguments,  k^ and  k^. 
Herterich  &  Hasselmann's  F(k.,t)  is  given  by  u>(kjC(k,t)/4gTr2  and  their  inter¬ 
action  coefficient  D  is  given  by 

2 

DfkgJCg.JOj )  =  -  (ujuj-j^Wj)  ^  T(k.,_k-j »kg jkg)  »  (3.4) 


for  resonating  quartets. 

The  above  identity,  Eq.  (3.4),  has  been  verified  numerically,  and  thus 
serves  as  a  mutual  check  of  the  rather  lengthy  algebra  involved  in  the  deriva¬ 
tion  of  both  models. 

Relation  to  the  nonlinear  Schrbdinqer  equation 

The  derivation  here  follows  the  lines  of  Zakharov  (1968),  who  showed 
that,  in  the  case  of  infinitely  deep  water,  the  cubic  Schrbdinger  equation 
is  a  particular  case  of  the  more  general  Zakharov  equation.  In  the  case  of 
finite  water  depth,  the  value  of  T(k_,k_1  .jcg ,k_3) ,  in  the  limit  when  .Jig *ii3 
tend  to  k  and  h  is  fixed  is  not  unique.  In  order  to  provide  a  better  grasp 
of  this  nonuniqueness,  we  include  here  an  outline  of  the  derivation  of  the 
finite  depth  nonlinear  Schrodinger  equation. 


-  15  - 


Restricting  the  analysis  to  narrow  spectra  around  k^  =  (kQ,  0),  we 
rewrite  all  wavenumbers  as  k.  =  +  ii»  *  (^,  X^)  and  |^|/kQ  «  1. 

Introducing  a  new  variable  A(^,t)  *  B(J<,t)e"^u(jL)'K!sp)]t  into  Eq.  (3.1) 
gives 


3A(*,t) 

1  — - [“(Jl)  “  w(k0)]A(^,t) 


w 

ll\  T(Jso*  t.  yti  »  JIq+412*  iio+^3^A  (^1 ) ACijjj ) A(^3 ) <5 (tjrhj/j -tj/g )  dh^di^cfy. 


(3.5) 


Equation  (3.2)  is  then  expanded  to  the  lowest  order  in  the  spectral  width 

i  ra)(k,)i  1/2  f  i[k  x1-u(krt)t]  7 

Mx.t)  *  2?  (tt-)  {e  |  A(i»^)e  +  c.c|  = 


Reja(x,t) 


iCk0x1-uj(ko)t3 


} 


(3.6) 


where  the  complex  wave  envelope  a(x,t)  is  the  inverse  Fourier  transform  of 
A(^,t),  multiplied  by  the  coefficient  (2w(kQ)/g)^2.  The  frequency  difference 
on  the  l.h.s.  of  Eq.  (3.5)  is  replaced  by  its  Taylor  series  expansion  up  to 
the  second  order  in  the  spectral  width 

c  c 1 

w(k)  -  wO^)  *  c  i|>  +  X2  +  -jf  $  +  0(  |ip|2)  (3.7) 

y  O  ~ 

where  cg  ■  3o)0/3k0>  c^  ■  32u)0/3k2  and  =  gkQth(koh).  Multiplying 

1  /2 

Eq.  (3.5)  by  (2wQ/g)  and  taking  its  inverse  Fourier  transform  yields 


c _  a2a  c'  a2a 

+  g  o  3  ,  _9.  3  A  _ 

2k  2  2  2 

c\  3Xg  c  3x^ 


fii  +  c  -ii-l 

\3t  g  sxj 

1  f2w  %l/2  ttt  i  (  + tb  rj  ~iiil  )  *  x 

!MT]  JjJ  T(^04^2+^3‘il ’V*i  *V^2’Vi3)A*(il dildi2d^3 

(3.8) 


Vt>. 


-%  A-V-.. 


» t  a. 


One-can  show  that  the  Taylor  series  expansion  ofT,  to  the  lowest  order  in  the 
spectral  width,  is  given  by 


-rtVVW  *'Ti  +  tii  *  Tm  +  T] 


(3.'9a) 


where 


T  °  -j.  [9o4  -10o2  +  9]  ;  o  s  th(k  h) 

1  3210*  0 

i*3  m  ...  \  2 , .  _  _  /■>  _2v ...  \2._/-,  2\2 


(3.9b) 


111  *SV 


k0  I  4cp  ( *  )  +4cpcq  ( 1  -°  )  ( )  4+9  ( 1  -ac )  *  lii  I  th  ( h  1&  -jM ) 


32ir  a  j«2  91^-^  lth(h|iJ-i1 1 )  -  • 

3  (3.9c) 

Y-  14-  -  12  +  l3o2  —  2o4]  (3.9d) 


1-4-  -  12  +  l3o2  -  2o4] 


[2c  +  cq(l-o2)2]2(V,1)2 


TIV - - B - * - f— 1 - o  ,  c  =  0,/kc 

1 6ir  .o  gth  ( h  |  jfeg- ^  | )  |  |  -  c2  ( ) 2  p 


(3.9e) 


Substituting  Eqs.  (3.9d,e)  into  Eq.  (3.8  )  gives 
fs.  ^  c_  .2.  c!  -2_ 


’(It*  V^)  +  “l|,|2a  +  cs  °z  u  (3‘10) 


where 


g2k4  kz  c 

al  * - 5  -  12  +  13a2  -  2o4)  ;  [2  ^  +  (1  -  o2)],  (3.11a,b) 

1 6u)q  o  o  g 


.  ?  i  (i9”Jb )  *2L 

“2  rr  (V^I )  A*(jh)A(jb)e 

i*-3)J— —  cz - r*i*  (3-lu> 

" 00  1^2“^  I )  "  g^-  (^2“^1^ 

Note'that  for  any  finite  depth  and  for  ■+•  0,  j  *  2,3  ,  the  values  of 
Tji,  Tjy  and  that  of  the  integrand  in  Eq.  (3.11c)  depend  on  angles  0j»  the 
"directions”  of  approach  to  the  limit,  where 


/.V  A.V.V.  %' ••a.- 


1 


-  17  - 


cose . 

J 


This  nonuniqueness  disappears  for  infinitely  deep  water. 

Regarding  the  integral  I,  Eq.  (3.11c),  one  can  show  that  I  = 
where  ^  is  a  solution  of  the  following  boundary  value  problem 


3X1  2=0 


for  z  <  0 


(3.12a) 


!!o  +  fi32*o  .  cq9°2  3 1 al 2 
31  9  aif  2“o  3xi  ’ 


at  2 


0 


3»o 

az 


0 


at  z  =  -h 


(3.12b) 

(3.12c) 


Thus_$0  appears  to  be  the  mean  flow  potential.  The  system  of  Eqs.  (3.10), 

(3.12)  was  obtained  by  Iusim  &  Stiassnie  (1982)  using  a  multiple-scale  approach. 
In  the  particular  case  where  the  water  depth  is  shallow  compared  to  the  group 
length  th(h|tj>j-i^  | )  can  be  replaced  by  an^  the  set  of  equations  given 

by  Davey  &  Stewartson  (1974)  is  recovered.  For  water  of  infinite  depth, 
Stiassnie  (1983)  extended  the  analysis  to  one  order  higher  in  the  spectral 
width,  and  rederived  Dysthe  (1979)  set  of  equations. 


This  section  deals  with  the  mathematical  formulation  of  one  of  the  , 
simplest  possible  non-trivial  nonlinear  interaction  problems  and  its  linearized 
(short  time)  solution.  The  smallest  number  of  wave  trains  required  to  enable 
significant  nonlinear  interaction  is  three  for  class  I  as  well  as  for  class  II 
interactions.  In  what  follows,  we  denote  these  3  waves  by  the  subscripts  a,  b, 
and  c.  For  anything  exciting  to  happen,  these  3  waves  have  to  form  a  nearly 
resonating  "quartet"  for  a  class  I  interaction,  and  a  nearly  resonating 
"quintet"  for  a  class  II  interaction,  see  Eqs.  ‘.(2.18)  and  (2.21),  respectively. 

To  form  a  "quartet",  or  a  "quintet",  out  of  three  waves,  one  can  "count" 
one  of  the  waves,  k  ,  twice  for  class  I  interactions,  and  three  times  for  class 

*  "a 

II  interactions. 

The  governing  equations  for  class  I  interactions  are  a  discretized  form 
of  Eq.  (3.1) 


dB  a  a  ^  ^  ifiyl 

’"dt*  "  <TaaJB/  +  2TabablBbl  +  2TacaclBc'  >Ba  +  2TaabcBaBbBce 


(4.1a) 


1  d T  ’  <2TI 


baba* 


I2  +  WBbi2  + 


2T«_  ■  |B  I  )B.  +  T.  B  B  e 
bcbc1  c1  '  b  bcaa  c  a 


*B2.'inIt 


(4.1b) 


dBr 

1  ^  -  (2T 


cacalBal2  +  2TcbcblBbl2  +  T< 


:lBcl  >Bc  +  TcbaaBbBae 


-iQjt 


(4.1c) 


where  8j  *  2u>a  -  -  wc. 

For  classll  3-wave  problems,  which  do  not  satisfy  Eq.  (2.18),  Eq.  (2.24) 
similarly  gives 


,!B,|2  +  2Tabab'Bbl2  +  2T 


acaclBcl2)B.  +  ^ialbc'OVc' 


1n„t 


(4.2a) 


-  19  - 


1  W-  ■  <2WBalZ  +  Tbbbb'Bbl2  +  2TbcbclB/>Bb  +  ub2cLBcBae’1!!nt 

'  (4.2b) 

1  ^  ’  '2Tc.calB.l2  +  2TcbcblBbl2  +  WUX  +  urtiaaBbBae”Snt 


(4.2c) 

(3) 

where  Qjj  »  3u>a-u>b-u>c  and  U^a'bc  is  assumed  to  be  synroetric  with  respect  to  b 
and  c. 

To  complete  the  mathematical  formulation  of  either  of  the  above  systems 
of  equations,  Eq.  (4.1)  or  Eq.  (4.2),  one  has  to  specify  the  following  initial 
conditions: 


Ba(0)  -  bfl,  Bb(0)  =  bb,  and  Bc(0)  =  bc  ; 


where  the  relation  between  Bi  and  the  actual  physical  amplitude  a-  is 


•j  ’  7  (i)  iBji 


One  can  assume  that  the  initial  amplitude  of  one  of  the  waves,  which  is 
called  "the  carrier"  and  denoted  by  the  subscript  a,  is  much  larger  than  the 
amplitudes  of  the  other  two  waves,  denoted  by  b  and  c,  to  be  called  "the 
disturbances":  |bb|,  |bj  «  |bfl|.  Only  linear  terms  in  the  disturbances  Bb, 

Bc  are  retained,  so  that  the  carrier  wave  remains  unaffected  in  this  short  time 

-1T,aa.lb«l2t 

analysis,  and  is  given  by  B„  *  bae  ;  b,,  is  assumed  to  be  real 

GO  O 

without  loss  of  generality. 

Class  I  Instabilities: 

The  wave  numbers  of  the  carrier  and  the  disturbances  are 


**  *=  0) 


«  kQ(l+p,  q),  =  k0(l-p;-q) 


(4.3a) 

(4.3b,c) 


■  . 


-K  -Le 


so  that  Eq.  (2.18a)  is  satisfied  identically.  The  linearized  version  of  Eqs. 
(4.1b,c)  is 


dBi  5  *  ,  “1nTt 

1  -dT  *  2Tbaba  baBb  +  TbcaaBcbae 


(4.4a) 


dB. 


2  *2 
/b_  +  T  . __Bub“e  1 


1  dt  =  2^caca  bauc  '  *cbaa“bua' 


(4.4b) 


where  flT  *  J1T  +  2T  b.  . 
l  I  aaaa  a 

'Assuming  a  solution  of  the  form 


Bb  •  v 


-i (0. Sfij+SjJt 


;  Bc  =  bce 


-i(0.5S2j-6j)t 


one  can  show  that  6j  must  be  given  by 

6 


i-  <Tbaba  -  Tcaca>  ba  1  D}/2 


(4.5a) 


where 


D,  -  [0.55,  -  (Tbaba  ♦  Tcaca)b|32  -  TbcaaTcbaabJ  (4.5b) 

Positive  values  of  Dj(p,q)  correspond  to  stability  regions  in  the  p,q  plane 
and  vice  versa.  The  curves  Dj  =  0  form  the  stability  boundaries,  and  the  point 
where  Dj  attains  its  minimum  is  called  the  most  unstable  mode.  The  value  of 
1  /2 

aI  *  (-Dj/gk0)  '  for  the  most  unstable  mode  is  called  the  maximum  growth-rate. 
Class  II  instabilities 

For  this  case,  the  carrier  wave  numbers  k^,  and  k^  are  still  given  by  Eqs. 
(4.3a,b)  but 


i£c  •  k0(2-P.  -q) 


(4.6) 


so  that  Eq.  (2.21a)  is  now  satisfied  identically. 
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The  linearized,  short-time,  version  of  Eqs.  (4.2b,c)  is- 
1  TT  ■  ^TbababaBb  +  “bcLa6^"’011' 


(4.7a) 


*e  2  (21  *  3  ~^II* 

i  “hT  *  2T^abaB,.  +  uih;aaBhbae 
dt  caca  a  c  cbaaa  b  a 


(4.7b) 


where  0,j  -  O,,  +  3Taaaal>a- 

Assuming,  again,  a  solution  of  the  form 


n  h  -i(°*5^n+,sii)t  -ito.SQjj-ajjjt 

Bb  bbe  »  Bc  "  bce 


one  finds  that 


SII  -  <Tbaba  '  Tcaca>ba  1  D1 If 


(4.8a) 


Dt,  »  [0.56,,  -  (T.  .  +  T  )b2]2  -  U<2>  U<2>  b 

UII  v'baba  'caca7  aJ  bcaaaucbaaa 


(4.8b) 


The  stability  boundary  and  maximum  growth  rate  for  class  II  interactions 
are  obtained  from  Eq.  (4.8b). 

Results 


In  Fig.  1,  we  show  the  class  I  and  class  II  instability  regions  (as  shaded 

zones)  for  kQh  *  2.  The  solid  lines  represent  the  calculated  results  and  the 

dashed  curves  are  those  of  McLean  (1982b, Figs.  2b  and  2c).  In  Fig.  la, 

k  a  *  0.195  (where  a  is  the  first  order  amplitude  of  the  carrier  in  Stokes' 

0  0  0 

expansion),  which  is  equivalent  to  (ka)m  *  0.2  (the  subscript  m  stands  for 
McLean).  As  a  conversion  formula,  we  used  the  following  expression 


(ka) 


24ch6(kQh)  +'3 
64sh6( kQh) 

given  by  Sbjelbreia  and  Hendrickson  (1961). 


m 


k«a-  + 
o  o 


(koao}  +  °(koao)' 


(4.9) 


Fig.  1  about  here  please 


In  Fig.  lb,  kQa0  «  0.326  (corresponding  to  (ka)m  *  0.35). 

The  locations  of  the  maximum  growth  rates  (pj,  qj),  (Pjj.  qj j ) ,  for  class  I 
and  class  II  instabilities,  respectively,  are  marked  by  x  for  our  results  and 
by  a  dot  for ‘McLean's  results,  and  their  numerical  values,  as  well  as  those  of 
the  maximum  growth  rates  Oj,  Ojj  are  given  in  the  figures.  The  overall  agree¬ 
ment  in  Fig.  la  (which  is  for  actual  amplitude,  which  is  471  of  the  theoretical 
maximum  amplitude,  Cokelet,  1977)  is  quite  satisfactory.  For  smaller  steep¬ 
nesses,  the  agreement  becomes  even  better.  On  the  other  hand,  for  very  steep 

waves  (in  Fig.  lb,  the  actual  amplitude  is  82X  of  the  theoretical  maximum 

/ 

amplitude),  the  agreement  is  less  impressive.  Nevertheless,  a  somewhat  better 
agreement  is  obtained  if  we  compare  McLean's  results  for  (ka)m  *  0.35  with  the 
artificially  amplified  value  aQk0  «  0.41,  see  Fig.  lc. 

A  similar  degree  of  agreement  was  obtained  for  several  other  water  depths, 
and  should  give  the  reader  some  indication  about  the  validity  of  the  present 

Fig.  2  about  here  p lease  J 

Figure  2,  which  is  quite  typical,  is  used  to  demonstrate  some  general  features 
as  well  as  clarify  some  of  the  terminology  which  is  used  later. 

One  can  see  that  a  certain  similarity  exists  between  class  I  and  class  II 
instability  regions.  Both  can  be  regarded  as  consisting  of  two  domains:  a 
wider  band  at  lower  values  of  p  and  usually  a  much  narrower  region  at  higher 
values  of  p.  The  first  region  will  be  referred  to  as  the  main  region,  and  the 
other  as  the  secondary  instability  region.  The  qualitative  difference  between 
class  I  and  class  II  instability  regions  is  that  for  class  I  the  two  domains 
are  usually  disconnected  while  in  the  case  of  class  II  they  are  bound  by  a  line 
of  Infinitesimal  thickness.  The  secondary  regions  sometimes  disappear  completely 
and  for  class  I,  the  instability  region  in  these  cases  terminates  at  some  q>  0 
(compare  with  Crawford,  et  al.,  1981,  for  infinite  water  depth). 

Figure  2  shows  the  three  wave-number  vectors  k^,  k^,  as  well  as  the 
location  of  four  points  of  local  maximum  growth  rates: 


;0(*? 


A.  Class  I  point  (pj,  qj)  with  local  maximum  growth  rate  Oj, 

B.  The  secondary.  Class  I  point  (p^,  q|)  with  local  maximum  growth  rate  Oj, 

C.  Class  II  point  (Pjj*  Pjj)  with  local  maximum  growth  rate  Ojj,  and 

D.  The  secondary  Class  II  point  (Pjj»  Qjj)  with  local  maximum  growth  rate  o*j. 

For  the  particular  data  of  Fig.  2  (kQh  *  0.35,  kQ^  =  0.04)  on  >  Oj  >  Ojj  >Oj. 

These  inequalities  are  by  no  means  general  as  will  be  shown  in  the  sequel. 
Nevertheless,  for  most  cases,  a**  >  a**. 

Fig.  3  about  here  please 

Figure  3  is  a  summary  of  the  results  for  Class  I  instabilities.  Figures 

3a  and  3b  give  the  values  of  Pj  and  q^,  respectively,  as  functions  of  the  water 

depth  and  wave  steepness.  The  depth  is  expressed  by  th(kQh)  (the  range  covered 

is  0.357  <  k  h  <  «);  and  the  wave  steepness  by  Cokelet’s  (1977)  e2,  denoted  here 
2  u  2 

by  e*  (the  range  0  <  <  0.7  is  covered). 

The  isolines  in  Figs.  3  and  4  were  drawn  using  interpolation  and  are  based 
on  about  forty  computed  data  points,  almost  equally  distributed  over  the  figure 
domain.  Figure  3c  is  a  plot  of  Oj  =  max(aj,  Oj)  isolines.  Note  that  for  the 
region  confined  by  the  broken  lines  a*  >  Oj  (sometimes  by  a  factor  of  three), 
whereas  the  opposite  is  true  in  the  outside  region.  For  the  case  where 
Oj  >  Oj,  p*  is  in  the  range  1.05t1.30  and  q j  =  0  which  implies  that  the  most 
unstable  mode  is  two  dimensional. 

Fig.  4  about  here  please 

The  results  for  Class  II  are  given  in  Fig.  4.  Note  for  this  case  p^j  is 
always  0.5.  Figure  4a  gives  the  values  of  qjj,  and  Ojj  is  shown  in  Fig.  4b. 

For  the  domain  above  the  dashed  line  in  Fig.  4b,  Ojj  >  o^»  which  indicates  that 
for  this  region  Class  II  instabilities  may  become  dominant. 

The  question  whether  the  disturbances  related  to  the  highest  value  of  o  will 
dominate  the  physical  process  remains  open,  and  awaits  additional  evidence.  The 
authors  hope  that  their  current  study  of  the  long  time  evolution  of  Class  I  and 
Class  II  instabilities  will  throw  some  light  on  this  and  on  other  relevant  aspects 
of  these  important  processes. 
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APPENDIX 


The  kernels  in  Eq.  (2.9) 
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The  interaction  coefficients  in  Eq.  (2.11) 
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Third  order: 


U'  1  /  *  U  _  U 

w0,l ,2,3  W1 ,2, -0,3  -0,1 ,2,3 


(A. 3a) 


0,1, 2, 3  ~  -0,-1, 2, 3  2, 3, -0,-1  2,-1 ,-0,3  -0, 2,-1 ,3 ~  -0,3, 2,-1  “ w3,-l ,2,-0 


W{3) 

W0, 1,2,3 
K(4) 

*0, 1,2,3 

where 


W0, 1,2,3 


Fourth  order: 

xO) 

A0,l  ,2,3, 
x(2) 

u,l ,2,3, 


v(3) 

*0,1 ,2,3, 


x(4) 

*0, 1,2,3 


(A. 3b) 

2W-0,-l,-2,3  "  W-0,3,-l,-2  +  W-l,-2,-0,3  "  2W-l,3,-0,-2  (A.3c) 
=  W0, 1,2,3  +  W1 ,2,0,3  (A*3d) 


■-fj  +  2li.,|th(|klh)- 

-  J  th(|^|h)th(|k,|h)  +  „>2+2  +  u^+3])  (A.3e) 


4  "  “0, 1,2, 3,4  +  B0,l, 2,3,4 


(A. 4a) 


4  =  a0, 4, 2, 3,-1  +  B0, 4, 2, 3,-1  +  a0,3,2,-l,4  +  B0, 3,2, -1,4  + 

+  “o,2,-l,3,4  "  B0,2,-l,3,4  ‘  °0,-l  ,2,3,4  '  B0,-l  ,2,3,4  (A.4b) 


4  =  a0,3,4,-lr2  +  B0, 3, 4, -1,-2  "  a0,-l,3,-2,4  “  60,-l,3,-2,4  " 

"  ot0,-l,4,3,-2  "  B0,-l  ,4, 3, -2  +  “o, 3, -2, -1,4  '  B0,3,-2,-l  ,4,  + 

+  °0, 4, -2, 3,-1  "  B0, 4, -2, 3,-1  "  °0,-l ,-2,3,4  +  B0,-l,-2,3,4 

(A. 4c) 

4  =  _a0,-l ,4, -2, -3  "  B0,-l ,4, -2, 3  +  “0,4, -2 ,-3,-1  "  B0,4,-2,-3,- 1  ~ 
"  a0,-l,-2,-3,4  +  B0,-l ,-2,-3, 4  “  “0,-1 ,-2,4,-3  +  60,-l ,-2,4,-3 


-Hi 


26  - 


t(5) 


X0,l ,2,3,4  *  "  “0,-1 ,-2,-3, -4  +  60,-l ,-2,-3, -4 


(A.4e) 


°o,i ,2,3,4 


1 

[<02“3“4'| 

32t? 

i2g“0^J 

1 

r“o“3“4] 

32tt3 

^guj-jujgj 

I  111  1 1 Jig 

O 

si. 

1 

tb(|k,|l.)  -  .l^n^thCL^IKHkjl]}  (A.4g) 


*The  kernels  of  Eq.  (2.15) 
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The  kernels  of  Eq.  (2.20) 
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laptions 


Fig.  1:  Bands  of  instability  for  kQh  «  2.  The  instability  boundaries 
are  given  by  the  solid  lines  and  the  points  of  maximum  growth 
rate  are  labeled  by  X.  McLean's  results  are  marked  by  the 
dashed  lines  and  by  • .  (a)  kQa0  ■  0.195,  (ka)ffl  *  0.2; 

(b)  kQao  «  0.326,  (ka)B  »  0.35;  (c)  k0a0  -  0.41,  (ka)m  =  0.35. 

Fig.  2:  Bands  of  instability  for  kQh  ■  0.35,  kQaQ  *  0.04  and  notation. 

Fig.  3:  Summary  of  results  for  Class  I  instabilities,  (a)  Isolines  of 

2  2 
Pj  (fOT  ec  *  0,  pj  *  0);  (b)  Isolines  of  qj  (for  either  c£  *  0 

or  th(kQh)  *  1,  qj  *  0) ;  (c)  Isolines  of  10 o™,  the  maximum 

2  jjj 

growth  rate  (for  ec  «  0,  Oj  *  0). 

Fig.  4:  Summary  of  results  for  Class  II  instabilities,  (a)  Isolines  of 
q  (b)  Isolines  of  lOo^. 
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